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Abstract 


A partial-wave analysis of A A elastic scattering data has been completed. 
This analysis covers an expanded energy range, from threshold to a laboratory 
kinetic energy of 2.5 GeV, in order to include recent elastic pp scattering data 
from the EDDA collaboration. The results of both single-energy and energy- 

dependent analyses are described. 
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I. INTRODUCTION 


This analysis of elastic nucleon-nucleon scattering data updates our previous analysis JTJ 
to 1.6 GeV in the laboratory kinetic energy. The present analysis extends to 2.5 GeV, which 
is the limit for clastic pp differential cross sections measured [|2[ by the EDDA collaboration 
using the cooler synchrotron at COSY. 

Measurements with a laboratory kinetic energy near 2 GeV are particularly interesting 
as they correspond to a center-of-mass energy (2.7 GeV) which has been suggested || for 
a clibaryon resonance Q. Near this energy, a sharp structure has been reported in the 
polarization observable A yy [||], and this was taken as support for such a resonance. A 
resonancclikc structure, at about the same energy, has also been reported in an analysis by 
Hoshizaki || . The authors of Ref. || have considered this possibility, but find no evidence 
for a resonant excursion in their cross sections. Polarization measurements expected from 
COSY and SATURNE II will certainly help to clarify this issue. 

The data base above 1.6 GeV is mainly comprised of cross section measurements, much 
of this coming from Ref. J^]. In Section II we describe the expanded database, noting the 
additions below 1.6 GeV as well as the new region from 1.6 GeV to 2.5 GeV. While the most 
significant changes are seen in our pp partial wave amplitudes, both pp and np data have 
been analyzed. 

In Section III, we briefly review the formalism used in our analyses. Here we present 
the updated amplitudes and make comparisons with our previous solution (SM94) [|l|. Fits 
with and without the new EDDA data are compared to show the influence of this particular 
measurement. Representative plots showing the agreement between our analysis (SM97) 
and cross section data have been generated to illustrate the quality of this fit. These results 
and the prospect for improvements are summarized in Section IV. 

II. THE DATABASE 

Our previous NN scattering analyses [|lj were based on 12838 pp and 10918 np data. In 
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Ref. 0 the pp analysis extended up to a laboratory kinetic energy of 1.6 GeV; the np analysis 
was truncated at 1.3 GeV. The present database [0 is considerably larger due both to an 
expanded energy range for the pp system and the addition of new data at lower energies. 
The total database is now about 20% larger than was used in our previous analysis [|T| . 

Below we list recent additions to our database. Some data sets which we collect are 
not used in the analyses, but are retained so that comparisons can be made. A complete 
description of the database and those data not included in our analyses is available from the 
authors [fjj. 

The new pp data have been produced mainly at COSY ||. From this source, we have 
added differential cross sections ranging from 540 MeV to 2520 MeV in the proton kinetic 
energy and from 35° to 90° in the cm scattering angle. In addition to this, about 60 
high quality polarized data (P, A xx , A yy , and A zx ) at 200 MeV were produced by the 
Indiana cooler ||. Another 35 high accuracy differential cross sections between 490 and 
790 MeV were recently published These measurements were made at LAMPF. We have 
added an excitation function of cross sections at 90° and between 0.3 and 0.4 MeV. These 
were measured at the Munster University low-energy machine |T0j. We have also added a 
measurement of A zz at 650 MeV produced by LAMPF [IT] but missed in the SAID database 

a. 

In constructing the data base extension from 1600 MeV to 2500 MeV, we reexamined 
a number of references in order to include higher energy data which had previously been 
neglected. This search netted additional data mainly from ANL (450 points) and Saclay 
(893 points). The complete set is listed Jl2| - |5l| in alphabetical order. 

The np database has not been increased significantly and, as a result, we did not extend 
our analysis of the 1 = 0 system. New np polarized data have been produced mainly by 


TRIUMF (101 points) ||52| , |53| , IUCF (33 points) [J54J, and LAMPF (49 points) |55| . The 
ANL-LAMPF-New Mexico University-Texas A & M University collaboration has hnalized 
its analysis of 311 high quality np polarized observables ( A xx , A ZZ1 A yy , and A zx ) between 
485 and 790 MeV and ranging from 25° to 180° ]56]]. These measurements were published 
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previously in Ref. |57|]. A few total cross sections in pure spin states between 4 and 16 MeV 
were produced by TUNL |5^| and Charles University at Prague []59| , |60| . Recently, the final 
LAMPF Any measurements between 480 and 790 MeV were also published |)T| . In addition, 
some new Ameasurements above 1190 MeV were made at JINR (Dubna) |62|]. Added 
unpolarized measurements include 15 np differential cross sections at 67 MeV from PSI 
and 6 differential cross sections at 14 MeV from Tubingen University |h|. A few missed 
differential cross sections at low energies from LAMPF |65| and at 1240 MeV from Berkeley 
H were also added. 

A few data sets were added to the data base but not included in the analysis. These 
include 82 missed np total cross section measurements between 4 and 231 MeV from LAMPF 
|67| . We excluded these data from the analysis in order to retain the same database (below 
350 MeV) as was used in the Nijmegen analysis This also applies to a new set of np 
differential cross sections at 162 MeV and at backward angles which were measured at the 
Svedberg Facility at Uppsala [|69 


III. PARTIAL-WAVE ANALYSIS 

Our first attempts to extend the range of the NN analysis used the parameterization 
scheme of Ref. [|l[ . These were unsuccessful. The problem was traced to the basis functions 
used to expand our K-matrix elements. Many of these become nearly degenerate as the 
kinetic energy of the incoming nucleon (T) increases to 2.5 GeV. As a result, a modified 
form was used in the present analysis. Apart from this difference, the formalism used here 
is identical to that used in Ref. [l|. The reader is directed to Refs. [7(].[71|] for more details. 
In the following we just outline the method used, in order to show how the modified basis 
functions fit into our parameterization scheme. 

For uncoupled partial waves ( 1 L ) 2 , 3 F 3 , ...), an S-matrix (S = SeSj ) is used. This product 
S-matrix is constructed from exchange ( Se ) and inelastic (S'/) pieces. Se is parameterized 
in terms of a K-matrix 
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S E = {l + iK E )/{l-iK E ), 


( 1 ) 


which in turn is expanded as 


K e = Born + a,Ai 


( 2 ) 


Here the Born term gives the singlc-pion exchange contribution and cq are free parameters. 
The expansion basis elements, An, are given by 


An = Fu 


T 


i —1 


where the function Fn, used as the expansion basis in our previous fits |l|, is given by 

'xq — x\ x l 


( 3 ) 


Fu(T) = 


4y(i 2 d / Tr\ - t\ F l/ 2 


Qi 


dx. 


( 4 ) 


MT Jo ^ ‘ V 1 — x J 1 — x 

Here M (/i) is the nucleon (pion) mass and xq = 1 + (4 /x 2 /MT). Qi is a Legendre function 
of the second kind. In Eq. (3), Tq is a parameter which was chosen to be 1 GeV. (The fit 
was not sensitive to this choice; fits using 0.5 GeV and 1.5 GeV were also attempted.) The 
basis function given in Eq. (4) was derived in Ref. JTT] . 

To ensure time-reversal invariance, the spin-coupled waves (for example, 3 P 2 — 3 F- 2 ) are 
parameterized as 


S(2 x 2) = S^SjS'J 2 , 


where again the matrix S E is expanded in terms of a K-matrix with the elements 


( 5 ) 


Km Born m H - ^ ^ Ai mi , 


( 6 ) 


the subscript (m = (+, 0, —)) labeling states with l m = (J + 1, J, J — 1). As in Ref. [jTJ , 
the matrix Sj is taken from a Chew-Mandelstam K-matrix coupling the NN channel to an 
appropriate N A state. This has been extensively described in Ref. 0- The simple modifi¬ 
cation of the basis elements, displayed in Eq. (3), provided the added flexibility required to 
extend our analysis to 2.5 GeV. 
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In Table I, we compare the energy-dependent and single-energy fits over the energy 
bins used in the single-energy analyses. Also listed are the number of parameters varied in 
each single-energy solution. A total of 144 parameters were varied in the energy-dependent 
analysis. 

Our single-energy and energy-dependent results for the isovector and isoscalar partial- 
wave amplitudes are displayed in Figs. 1 and 2. Here we also compare with our previous fit 
(SM94). In some cases the changes are quite large. This is particularly true near the upper 
energy limit of SM94, and for the smaller partial waves. The effect of these changes can be 
clearly seen in Fig. 3, where we show how well the new EDDA data @ are reproduced by 
both SM94 and SM97. The influence of this experiment is most pronounced in the forward 
direction. 

In general, we find little structure over the higher energy region. This reflects the smooth, 
and rather flat, total and reaction cross sections between 1.5 GeV and 2.5 GeV. Our fit to 
these quantities is displayed in Fig. 4. Note that the reaction cross sections were excluded 
from our fit. This verifies that the set of total, total elastic (deduced from differential cross 
sections), and reaction cross sections are self-consistent. 

The present analysis actually gives an improved fit to the data below 1.6 GeV. This 
is due to the altered basis set, found necessary to fit the higher energy data. Numerical 
comparisons are given in Table II. Here we see that the COSY data || comprise a large 
fraction of the total set above 1.6 GeV. The results of analyses with (SM97) and without 
(NM97) this data set show how influential these measurements have been in determining 
the amplitudes. (The fits SM97 and NM97 used identical parameterization schemes. Only 
the data base was changed.) The COSY data contribute a y 2 /datum of 1.07 when included 
in the fit. This jumps to 5.6 when we attempt a prediction based on the remaining data. 

IV. CONCLUSIONS AND FUTURE PROSPECTS 

We have extended our pp partial-wave analyses nearly 1 GeV beyond the limit quoted 
in our previously published results |l|]. The present range has been selected to include all 
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of the recent elastic pp cross section data measured by the ED DA group 0. We found 
that it was possible to simultaneously fit the pp total cross section data, in particular the 
precise data of Ref. |23|], along with differential cross sections from the EDDA collaboration 
[2|. The resulting reaction cross sections, which were not included in the fit, are quite well 
reproduced. The predicted reaction cross sections are consistent with the results of Ref. |7 2 


at lower energies, but deviate from these and follow the results of Ref. |73| above 1 GeV. 

While we find that the partial-wave amplitudes above 1.6 GeV are smooth and struc¬ 
tureless, reflecting the behavior seen in the total and elastic cross section data, we have also 
considered the effect of more localized structures reported in polarization measurements 
[jp. We can add resonancelike structures in individual partial-waves to see their effect on 
any observable. This will be utilized as more polarization data become available. 

As the high energy region was constrained mainly by cross section data, the present solu¬ 
tion should be considered as a guide to the expected amplitudes. The EDDA collaboration 
is planning to measure P, A yy , A xx , and A xz in the near future. This will be crucial to any 
future analyses. 

Further data is also expected from a number of other labs. About 2000 polarized pp 
measurements are expected above 1000 MeV [75] as the nucleon-nucleon program at SAT- 
URNE 11 is completed. While not included in the present fit, preliminary data j75| from 
SATURNE If is in reasonable agreement with our predictions. A representative fit to P 
data, at 2.16 GeV, is given in Fig. 5. 

A similar number of polarized quantities from np elastic scattering are expected (between 
250 and 560 MeV ]]7(J) from PSI. The Freiburg University group is also planning to publish 
np measurements which were done at PSI at the beginning of the 1980’s. These data 
range from 200 to 580 MeV and from 77° to 179° |]77[]. Final np differential cross sections 
between 73° and 179° measured at Uppsalla ][78| are expected to replace data at 96 MeV 
79j] and 162 MeV ][d 9|]. IUCF is also measuring np differential cross sections in the backward 


direction at about 200 MeV to solve a shape problem in the angular distribution [|5D|. Other 
np sources include an extension of A ul measurements |R| at JINR [^], TRIUMF analyzing 
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power measurements at 350 MeV |82[ , and TUNL measurements |83|] of the P parameter and 
A a L at 7 and 15 MeV. We will continue to update our energy-dependent and single-energy 
solutions as the new measurements become available. 


Finally we note that by extending our analysis to 2.5 GeV, we may be bridging the gap 
between the low- and high-energy regions. This is suggested if we plot da/dt versus s, as is 
shown in Fig. 6. The result expected from dimensional counting at high-energy and fixed 
cm angle is 


da 

dt 


1 

pvW 



(7) 


where N is the minimum number of fundamental constituents (quarks). While a slightly 
extended energy range would be more definitive, our results do appear to be consistent with 
this limit. 
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FIGURE CAPTIONS 


Figure 1. Isovector partial-wave amplitudes from 0 to 2.4 GeV in the proton kinetic 
energy. Solid curves give the amplitudes corresponding to the SM97 solution. 
The real (imaginary) parts of the single-energy solutions are plotted as triangles 
(squares). For comparison, the previous solution SM94 JI| is plotted with (+) 
marks. The (x) marks give Im T-T 2 -T^ from SM97, where T ^ is the spin-flip 
amplitude. All amplitudes are dimensionless. 


Figure 2. Isoscalar partial-wave amplitudes from 0 to 1.2 GeV. Notation as in Fig. 1. 


Figure 3. Comparison between SM97 (solid curve) and differential cross sections at 
(a) 9* = 45° ± 1° and (b) 9* = 90° ± 1°. Recent COSY measurements [|| are 
plotted as filled circles. Other data from the SAID data base [0 are plotted as 
crosses. Our previous solution (SM94) is plotted to 1.6 GeV (dot-dashed line). 


Figure 4. Total cross section comparisons, (a) The solid (dashed) curves give the pre¬ 
dictions of solution SM97 for the total (total clastic) cross section. Experimental 
points are from the SAID data base [[71]; filled circles are from Ref. ||23||. (b) 


The solid curve gives the total reaction cross section of SM97. Filled circles are 


estimates from Ref. |72j. Filled triangles are estimates from Ref. 


Figure 5. Angular dependence of recent SATURNE II analyzing power (P) data [75 


This measurement, at 2.16 GeV, was not included the SM97 analysis. The solid 
line gives the SM97 prediction. The dashed lines are generated from a single- 
energy solution and its associated error estimate. 


Figure 6. da/dtt plotted as a function of s at 9* = 90°. The SM97 solution is plotted 
as a solid curve. The dash-dotted line gives da/dt ~ s -10 . The plotted data are 


from Ref. [86 
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Table I. Comparison of the single-energy (SES) and energy-dependent (SM97) fits to pp 
and rip data. Values of y 2 are given for the SES and SM97 fits (evaluated over the same 


energy bins). Also listed is the number of parameters varied in each single-energy solution. 


Energy Range (MeV) 

X 2 SES(SM97)/pp data 

X 2 SES(SM97)/np data 

Parameters 

4-6 

22(39)/28 

50(66)/53 

6 

7-12 

84(132)/88 

221(309)/87 

6 

11-19 

17(49)/27 

191(445)/236 

8 

19-30 

123(275)/114 

263(286)/295 

8 

32-67 

294(375)/224 

667(754)/485 

10 

60-90 

55(64)/72 

457(595)/329 

10 

80-120 

161(185)/154 

419(487)/353 

10 

125-174 

301(310)/287 

328(367)/272 

11 

175-225 

249(354)/212 

715(766)/499 

13 

225-270 

66(91)/64 

243(270)/236 

13 

276-325 

274(309)/256 

571(655)/518 

17 

325-375 

297(320)/246 

421(474)/353 

17 

375-425 

555(601)/436 

753(843)/549 

17 

425-475 

902(1004)/665 

775(799)/629 

18 

475-525 

1322(1484)/1081 

1252(1419)/787 

30 

525-575 

861(972)/754 

549(584)/432 

31 

575-625 

1032(1154)/760 

422(491)/367 

36 

625-675 

891(863)/754 

1263(1563)/875 

36 

675-725 

838(882)/777 

403(473)/386 

37 
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Table I. (continued) 


Energy Range (MeV) 

X 2 SES(SM97)/pp data 

X 2 SES(SM97)/np data 

Parameters 

725-775 

990(1195)/827 

512(558)/374 

37 

775-824 

1583(1754)/1170 

1518(1845)/944 

38 

827-874 

1195(1358)/939 

386(497)/366 

39 

876-924 

341(412)/389 

753(920)/628 

41 

926-974 

790(945)/679 

354(498)/352 

43 

976-1020 

931(1131)/708 

300(441)/331 

43 

1078-1125 

528( 689)/413 

427(671)/326 

45 

1261-1299 

680(972)/507 

— 

29 

1481-1521 

139(266)/149 

— 

29 

1590-1656 

472(655)/409 

— 

31 

1685-1724 

185(293)/118 

— 

31 

1778-1818 

404(628)/347 

— 

31 

1929-1968 

218(271)/168 

— 

31 

2065-2104 

673(1241)/431 

— 

31 

2176-2224 

1005(1325)/377 

— 

31 

2330-2470 

803(1257)/458 

— 

31 
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Table II. Comparison of present and previous solutions. Dataset A was used in the SM94 
analysis [jTJ . Dataset B contains all data (apart from the EDDA data 0) used in generating 
solution SM97. See the text for details regarding the SM97 and NM97 fits. 


PWA 

Data 

X 2 /pp data 

X 2 /np data 



(0-1600 MeV) 

(0-1300 MeV) 

SM94 0 

(dataset A) 

22375/12838 

17516/10918 

SM94 0 

(dataset B) 

22390/12889 

18480/10843 

SM97 

(dataset B) 

20910/12889 

17400/10843 



(0-2520 MeV) 

(0-2000 MeV) 

SM97 

(dataset B) 

26460/14873 

17440/10854 

SM97 

(EDDA dataset [||) 

2278/2121 

— 

NM97 

(dataset B) 

25240/14873 

17280/10854 

NM97 

(EDDA dataset @) 

11964/2121 

— 
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